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■ We derive the first six coefficients of the heat kernel expansion for the electro- 

magnetic field in a cavity by relating it to the expansion for the Laplace operator 
acting on forms. As an application we verify that the electromagnetic Casimir 
energy is finite. 
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1 Introduction 



X 



The modes of an electromagnetic field in a cavity, taken together with their unphysical, 
longitudinal counterparts, can be mapped onto the eigenstates of the Laplacian acting on 



the de Rham complex of a 3-manifold with boundary. The electric and magnetic fields are 
thereby associated to forms of degree p — 1 and p = 2 respectively. In this correspondence 
transverse modes are associated with coexact, resp. exact forms, which permits to further 
map longitudinal modes to forms of degree p = and p = 3. We will use this observation, 
which is explained in detail in Sect. |21 below, to compute the first six coefficients of the 
heat kernel expansion for the electromagnetic field in a cavity. The result is used to show 
in a simple way that the Casimir energy in an arbitrary cavity with smooth boundaries is 
finite, a conclusion which has been reached previously [B]. In an appendix the derivation 
of the numerical coefficients of the expansion is presented. 

We shall present a Hilbert space formulation of the classical Maxwell equations in a 
cavity Q C M 3 . In a preliminary Hilbert space L 2 (f2,IR 3 ) we define the dense subspaces 

TZ = {V G L 2 (fi,M 3 ) | rotV G L 2 (fi,M 3 )} , 
U = {V G TZ I (U,rotV) = (rotU,V), VU G 11} 
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and the (closed) operator 

R = rot with domain T>(R) = TZq . 

Its adjoint is then given as R* = rot with T>(R*) = 1Z. We remark that R, resp. R*, is 
also the closure of rot defined on smooth vector fields V with boundary condition Vn = 
on the smooth boundary dQ, resp. without boundary conditions. This is what is meant 
when we later simply say that a differential operator is defined with (or without) a certain 
boundary condition. 



The subspace 

H = {V G L 2 (ft,M 3 ) | divV = 0} (1) 

and its orthogonal complement in L 2 (Q, IR 3 ) are preserved by R and, therefore, by R*. We 
will thus view them as operators on the physical Hilbert space TC. The Maxwell equations 
with boundary condition Eii = on the ideally conducting shell dQ can now be written 

as 

' E 



d 


(» 




-iR 


iR* 




M[ b ) (2) 



with 

m = i , i . m <>n n + n 



cf. ^21- Since no boundary condition has been imposed on B, we have M(0, B) = for 
all B = V"0 with ij) harmonic, and hence 

dim Ker M = oo . (3) 

We shall compute the heat kernel trace 



-twt 



where ' means that the contributions of zero-modes, i.e., of eigenvalues Uk = of M, have 
been omitted. This is necessary in view of (JBJ), but a more physical justification, tied 
to the application to the Casimir effect to be discussed later, is that zero-modes are not 
subject to quantization. 



The square of M is 

M 2 = (R*R o W-a e o \ 

V RR* J \ -A B J ' 1 ' 

where Ae, resp. Ab, is the Laplacian on 7i with boundary conditions 

E||=0, resp. (rotB)||=0. (5) 

The operators RR* and R*R have the same spectrum, including multiplicity, except for 
zero-modes. Incidentally, we note that eigenfunctions (E, B) corresponding to Uk ^ 
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satisfy B = — ia;^ rotE and hence, by Stokes' theorem, the boundary condition Bj_ = 0, 
which we did not impose, but which is usually also associated with ideally conducting 
shells. Since <9 f 2 + M 2 = (id t — M)(—id t — M), each pair of non-zero eigenvalues of R*R 
and RR* corresponds to a single oscillator mode for (j2j). We will thus discuss the heat 
kernel asymptotics for 

I TV (e~* M2 )- J Tr « e * AE ( 6 ) 
2 i] W e j -\ T 4e tA B (7) 

oo 

-^a n t^, (HO). (8) 

n=0 

The coefficients a n are known, see e.g. for general operators of Laplace type. The 
direct application of such results is prevented by the divergence constraint in H, see (JTJ). 
In the next section we indicate how to remove it. First however we present the main result. 

Let 

L a6 = (V ea e t ,n), (a, b =1,2), 

be the second fundamental form on the boundary dfl with inward normal n and local 
orthonormal frame {ei,e2,n}. We denote by \Q\ the volume of Q and set 



f[dQ] = [ f(y)dy 
Jan 



where dy is the (induced) Euclidean surface element on dfl. The corresponding Laplacian 
on dfl is denoted by V 2 . 

Theorem 1 Let Q C 1R 3 be a compact, connected domain with smooth boundary dQ 
consisting of n components of genera gi, g 2 , . . . , g n . Then 

a Q = 2(4tt)-#|0| , 
a,\ — , 

a 2 = -~(47r)-§(trL)[c>Q] , 

1 1 n 

a 3 = -(47r)- 1 (3(trL) 2 -4detL)[^]--^(l + ^) + l, (9) 

i=l 

16 3 

a 4 = (47r)-5(2(trL) 3 -9trL-detL)[<9fi] , 

a 5 = — I (47r)- 1 (2295(trL) 4 - 12440(tr L) 2 det L + 

122880 

+ 13424(detL) 2 + 1200trL • V 2 trL)[d^] . 

We will give two partially independent proofs, based on (JHJ), resp. (0). Their agreement 
is related to the index theorem, as it may be seen from (£Q). A further, partial check of 
these coefficients has been made on the basis of general cylindrical domains and of the 
sphere, where a separation into TE and TM modes is possible. 
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The coefficient a was computed in [T3] (except for the factor 2 replaced by 3, as the 
divergence condition ((TJ) was ignored), a 1; a 2 in pQ. The coefficient a 3 is closely related to 
a result of [3 , as discussed in Sect. H3 



2 Proofs 

We consider the space of (square integrable) forms, A(f2) = 0™ =o A P (Q), on the manifold 
Q with boundary, together with the exterior derivative d p+ \ : A P (Q) — > A p+ i(Q) defined 
with relative boundary condition ([Hj, Sect. 2.7.1) 

as a form uj\qq G A p (dQ). For later use we recall that by the de Rahm theorem for 
manifolds with boundary (|J] or [TI], Thm. 2.7.3) we have 

H p r (Sl) = # n _ p (Q) S , (10) 

where H^(Q) = Ker d p +\/Im d p is the p-th relative cohomology group, H P (Q) is the p-th 
homology group, and H p (Q,dfl) is the p-th relative homology group, i.e., the homology 
based on chains mod dQ. 

We shall henceforth restrict to Q C M 3 as in Theorem ^ Using either homology 
the dimension of H^(Q) is seen to be 






(p = o) , 


n — 1 


(p = 1) , 


i=l 


(P = 2) , 


1 


(p = 3) . 



(11) 



These are also the dimensions of the spaces of harmonic p-forms. 
The space A(Q) = ©p =0 A P (Q) may be identified as 

A(O) = L 2 (ft) © L 2 (n, R 3 ) © L 2 (Vl, R 3 ) © L 2 (ft) 9 (0, E, B, ^) , 
where d : A(f2) — > A(f2) acts as 

d : L 2 (ft) — ► L 2 (ft,R 3 ) — ► L 2 (ft,R 3 ) — ► L 2 (Q) — >0 

grad rot div 

with boundary conditions = 0, Em = 0, = on dQ. Then 

d* : < — L 2 {Q) < — L 2 {Q,R 3 ) < — L 2 (fi,R 3 ) < — L 2 {Q) 

—div rot —grad 
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without any boundary conditions. The Laplace-Beltrami operator on forms, 

3 

-A = 0(-A p ) = dd* + d*d , 

p=0 

is seen to correspond to the Euclidean Laplacian with boundary conditions 

= (p = 0) , 

E||=0, divE = (P=l), 

B ± = 0, (rotB)||=0 (P=2), 

(grad^)j_ = (p = 3) . 

Each of the four problems admits a heat kernel expansion, 



(12) 



OO 

n — 3 



Tr M n)e^ = $>i^, (13) 

n=0 

whose coefficients have been computed (n — 0, ... ,3) [3] or can be computed using existing 
results (n — 4,5) [3]. To this end we note that the boundary conditions for p = 1, 2 can 
be formulated equivalently as 

<9E 

E,|=0, — ±-(trL)E ± = (p = l), 

0n (M) 



<9B 



B ± = 0, -^J1-LB|| = (p = 2). 

First approach. We will compute (JUJ). We observe that — Ae is just the restriction of 
— Ai to its invariant subspace 

H = {Ee L 2 (n,R 3 ) | divE = 0} = Ker d\ . 

Hence 



TV' P *A E _ rp a tAi _ ijl/ 

ir w e — ir L 2( n M 3) e ir w± e 



H ± t 

where the orthogonal complement of in L 2 (f2, R 3 ) is 

^ = Ran ck = Ran di = { V0 G L 2 (fi, R 3 ) | = on dQ] , 

(Ran d is closed by the Hodge decomposition, see e.g. [HUH))- By dA = Ad, the operators 
(—Ax) \ n ± and — A have the same spectrum (in fact V0 = implies = by the 
boundary condition). Thus, using also (JTTJ), we find 

TV J A E — <Vr' J A 1 TV' o* A 

Lr H e — ir L 2 (f2R 3) e - ir L2(n) e 

= Tr L 2 (CiM3) e' Al - Tr L 2 (c) e* A ° - (n - 1) , 

i.e., 

flfc = 4 1} - 4° } , (k ^ 3) , 
(i) (o) , -i 

a 3 — a 3 a 3 W + 1 . 
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These relations, together with the values of computed in the Appendix, yield the 
values of the coefficients stated in the Theorem d In particular, we will obtain 

I 



3 



— (4vr)- 1 (3(trL) 2 + 28detL)[9fi] . 
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This matches the stated value of a 3 because of 

- n 1 n 

i=l i=l 

and of the Gauss-Bonnet theorem, 

1 n 1 

-^(l-^) = 2( 47r ) _1 ( detL )[ a °]- ( 15 ) 

i=l 

Second approach. We now compute ((Zj). As has been noted in the Introduction, 
eigenmodes of — Ab, except for zero-modes, satisfy the boundary condition = 0, and 
are thus eigenmodes of — A 2 belonging to its invariant subspace H, cf. (0 fT2|) . The 
converse is obvious. We conclude that 

li n e — 1^2(^3) e - lr n ±e 

Since 

H = {B G L 2 (fi,M 3 ) | divB = 0} = Ker rf 3 , 

we have 

ft x = Ran d* 3 = Ran d£ = {-V^ G L 2 (ft,R 3 ) | ^ G L 2 (fi)} . 

Using <i*A = Ad*, we see that (— A 2 ) and —A3 have the same spectrum, except for 
a single zero-mode (in fact, — = implies tj) = const ). We thus find, using (fTTj) . 



rp I tA B _ rpl tA 2 rp / tA 3 

lr w e - ir L 2 (f2iR 3) e - ir i2(n) e 

rt 

= Tr L 2 (Q R3) e tA2 - Tr L 2 (n) e* Aa - Q^ft - l) 



i=i 



i.e., 



a fc = 4 2) - 4 3) > ( k ^ 3 ) > 

n 

(2) (3) 11 

03 = a 3 ~ a 3 - 9* + 1 • 

1=1 

From these relations and from the results of the Appendix we again recover Theorem ^ 
In particular, 

4 2) -4 3) = ^(47r)- 1 (3(trL) 2 - 36 det L) [<9ft] 
leads to the claim for 03, because of 

n 1 n 1 n 

i=l i=l i=l 

and of (JED. 
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3 Application to the Casimir effect 



For the purpose of this discussion we simply define the Casimir energy by the mode sum- 
mation method, see e.g. In particular, we do not address the issue jU] of whether it 
is the most appropriate physically. We shall however observe that the Casimir energy is 
finite - a conclusion obtained in [3], but questioned in [TU] . 

Consider the cavity Q C 1R 3 enclosed in a large ball Qq. As usual we compare the 
vacuum energy of the electromagnetic field in the domains Q U (Qq \ Q) with that of the 
reference domain Q . Each eigenmode of either domain contributes a zero-point energy 
Uk/2, resp. uj®/2. As a regulator for the eigenfrequencies Uk = A^ 2 , we choose e _7Afc , 
(7 > 0). The corresponding definition of the Casimir energy is 



E c =\\sm lim(j]A|e-^ - 5>2)*e 

Z Slo^oo 7J.O \ * — ' L — ' 
\ k k 



We shall prove that the limit 7 j is finite. It will also be clear that the subsequent limit 
Q — > 00 exists, though we shall not make the effort to prove that (see however e.g. jH], 
Section 12.7 for the necessary tools). Using 



id 
dt t"2— e 
dt 



-tXk 



and (JEJ) we find for the regularized sum of the eigenfrequencies 



£4 



n 



2V^ an Jo 



k n=0 

as 7 I 0. Here 5 > is arbitrary, but fixed, and 



dt £~2(t + 7) 2 



means up to terms 0(1). Using 



dt t"2(£ + 7 )V 



4 — 2 



f7~ 

27- 1 
_i 

7T7 2 

-log 7 



(n = 0) 
(n = 1) 
(n = 2) 
(n = 3) 
(n = 4) 



we find 



E4 



2 e - 7 A fc 



a 7" 



7T _3 

017 2 



1 



a 2 7 1 + ■ a 3 7 2 + — ^a 4 log 7 

77 2V7T 



Hence a finite Casimir energy requires (cf. |7j) that CLQy (3^, 0,2; ^4 

(but not necessarily 

03!) agree for Q U (Qq \ ^) and for the reference domain Qo- This is indeed so for 

3 

a = 2(47r)~2 \Q \ and for a% = 0, but also for 02, 0-4 as the contribution from the two 
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sides of dQ cancel. The same conclusion is obtained if the regulator e 7Afc is replaced by 
e -(7A fc ) 1/2 ( see |7 ; E q. (27)): 

1 12 24 3 2 i 1 

y~] A|e~ (7Afc)1 « ^a 7~ 2 + 4ai7~2 + — a 2 7 _1 + ■ a37~" + -^a 4 log7 . 

y7T y'Tl y'TT 

Since no renormalization is necessary, the value of Eq agrees with that obtained by means 
of the zeta function. 

In the rest of this section we compare our results with those of |U El To the extent 
the comparison is done we will find agreement. An important tool there is the mode 
generating function, Eq. (4.5) in [2], 



2 \-A E -k 2 -A B -fc 2 
= j Tr'((-A E - k 2 )- 1 + ("Ab - k 2 )- 1 ) , (keC\R) 



(16) 



where "=" means equality "within addition of some polynomial in k 2 " . Since the resol- 
vents in (|16j) are not trace class, but their squares are, we first consider that replacement. 
Using (A + [i)~ 2 = J °° dt t e~*( A+/i ) we obtain, as // — > oo, 



"1 OC it OO 

iTr'((-A E + Ai r 2 + (-AB + /i)- 2 ) =J2 a n / dt-t^e-* = 5>(^)a n/ ^ 

n=0 ^° n=0 

with coefficients a n given in Theorem ^ Integrating w.r.t. /i we find 

1 oo 

- Tr' ((-A E + /i)- 1 + (- A B + /x)" 1 ) = £ r^K/i"^ - a x log// 

n=0 

and hence, with /i 1 / 2 = — ik, 

$(k) = 2^a ik 3 - y/naik 2 ln(-A; 2 ) + \^Ka 2 k - a 3 + 0(kr x ) . 

Upon insertion of the mentioned values for do, . . . ,0,3 this agrees with Eq. (4.40) in 
except for 03 which is there replaced by its local part, see ©, 

~a 3 = ^(47r)^(3(trL) 2 -4detL)[^] = ^ J da(^{K 2 + k 2 2 ) - Ki k 2 

where Ki, k 2 are the principal curvatures. Note however that this discrepancy is implicit 
in the definition of "=". It is resolved in [3] by first considering 5$>(k), i.e., the difference 
of the mode generating functions corresponding to the configurations Q U (SIq \ fl) and 
Q . Thus 

6Q(k) = -2~a^ + 0{k- 1 ) , 
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since the contributions to a , cancel, and those to 03 double the value. Not ambiguous 
then is "the number of additional modes of finite frequency created by introducing the 
conducting surface <9f2" : 

C = 0(0+) - 0(oo) , 

where ip(y) = 5&(iy). For a connected boundary dQ of genus g the value of '0(0+) has 
been established as ip(0+) = —g (see Eq. (5.8)), resulting in 

C = 2a 3 -g. (17) 

This result agrees with Theorem ^ the non-local terms in take the values — 1(<? — 
1), — ~g, \ for fl, Q Q \ Q and flo respectively. Thus, 

8a 3 = 2a 3 - g , 

in agreement with (fTTj) . 



A Appendix 

In this appendix we compute the heat kernel coefficients in ()13|) for p = 0, . . . , 3 and 
n = 0, . . . , 5 on the basis of Theorems 1 and 4 in j^]. We use the same notation, together 
with P = n ® n denoting the normal projection at the boundary. The vector bundle 
is 1/ = n x R for p = 0,3, resp. V = TQ for p = 1,2, equipped with the Euclidean 
connection. The decompositions of V\q^i = Vn ® Vp 3 (4> N ,4> D ) (with projections n + , 
resp. Il_) and boundary conditions (p N n + S(fi N = 0, resp. 4> D = 0, are specified as follows, 
cf. (HI and 0: 

p = : 
p — 1 : 
p = 2 : 
p = 3 : 



n + 


= 0, 






n 


= 1, 






n+ 


= p, 


s 


— ~L aa P 1 


n 


= 1 - 


p, 




n+ 


= 1 - 


p, 


S = -L , 


n 


= p, 






n+ 


= 1, 


S -- 


= 0, 


n 


= . 







;is) 
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The result is 

a? = \(A7r)-^\dQ\, 

4 p) = i(4vr)-l4 p) (trL)[^], 

4 P) = ^ i (47r)- 1 (4?(trL) 2 + C3 ?(detL))[^] , 

a? = -i-(47r)-i (cg(tr L) 3 + eg tr L ■ det L) [02] , 

4 P) = ^^( 47T )^ 1 ( c 5i ) ( trL ) 4 + 4?(trL) 2 detL + c^detL) 2 + cgtrL ■ V 2 trL) [Ol] 
with coefficients given by 





p = 


p — 1 


p = 2 


p = 3 


JP) 
(p) 

c l 

JP) 
c 2 


1 

-1 
1 


3 

-1 
-3 


3 
1 

-3 


1 
1 
1 


JP) 
c 31 


3 


21 


33 


15 


>) 
°32 


-20 


148 


-220 


-4 


>) 
t-41 


4 


36 


60 


28 


Jp) 
c 42 


-18 


-162 


-186 


-42 


Jp) 
c 51 


555 


5145 


8625 


4035 


>) 
c 52 


-2840 


-27720 


-35720 


-10840 


Jp) 
c 53 


2224 


29072 


29712 


2864 


>) 
c 54 


120 


2520 


4680 


2280 



These values imply Theorem ^ as explained in its proof. 

The computation of the table is based on the general result of jS], which has been 
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applied to (|T5j) using the following identities: 



Tr(P a P 6 ) = 2(L 2 ) ab , 
Tr(P :a P a P :b P :b ) = (L 4 ) aa + (L 2 ) aa (L 2 ) bb , 
Tr(P :o P 6 P a P 6 ) = 2(L A ) aa , 

T , a/tA\ , a/t2\ t t1\ 

>bb 



Tr{P :aa P :bb ) = 2L ac:a L bc:b + 4(L 4 ) aa + 4(L 2 ) aa (L 2 " 



^(P-.abP-.ab) — 2L a b :c L a b :c + 6(L 4 ) aa + 2(L 2 ) aa (L 2 ) bb . 

They can be derived by using V ea n = —L ab e b , so that 

P- a = -L ac (e c <g> n + n <g> e c ) , 

and by assuming without loss that V ea e b has no component parallel to T p dQ at the point 
p of evaluation, i.e., V ea e b = L ab n. Then 

P :ab = -L ac . b (e c ® n + n <g> e c ) - 2(L 2 ) ab P + (L ac L bd + L ad L bc )e c <g> e d , 

from which the above traces follow. In turn they allow the computation of similar traces 
with P replaced by x — n+ — n_, i.e., by x — ±(2P — 1) in the cases p = 1,2. In these 
two cases we also have 

Tr S- a = —L bb:a , 
Tr S :ab = L cc:ab , 

and, moreover, for p — 1, 

:a) L bb - a L cc - a 

+ 2L bb L cc (L ) 
Tr(P :a S' : 6) = —2(L 2 ) ab L cc , 

L bb :aL C c:a ~\~ L bb L cc (L ) a a j 



resp. for p = 2, 

T^{S : aS- a ) = Lab:cL a b:c + 2(L 4 ) aa , 

Tr(P :a 5 :a ) = 2(L 3 ) aa , 
Tr(PS :a S :a ) = (L 4 ) aa . 

Furthermore, traces of L fc , (/c > 2), were reduced to tr L, det L by means of L 2 — (tr L)L + 
det L = 0. Finally, we used the Codazzi equation, L ab -. c = L ac:b , as well as 

L a b:ca — L ab ;ac = L aa (L 2 ) bc — (L 2 ) aa L bc , 

which follows from the Gauss equation. 
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